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SUMMARY 


A study has been made of nitric oxide formation in a laminar CO-air diffu- 
sion flame over a pressure range from 1 to 50 atm. The carbon monoxide (CO) 
issued from a 3. 06-mm-diameter port coaxially into a coflowing stream of air 
confined within a 20. 5-mm-diameter chimney. Nitric oxide concentrations from 
the flame were measured at two carbon monoxide (fuel) flow rates: 73 standard 

cm 3 /min (seem) and 146 seem. 

As the pressure was increased above i atm, the flames changed shape from 
wide and convex at the lowest pressures to slender and concave at 50 atm, 
with the most noticeable change occurring below about 10 atm. Flame height 
decreased as the pressure was increased from 1 atm to about 1 0 atm, but height 
was independent of pressure at higher pressures. The regimes of stable burning 
decreased with increasing pressure up to 50 atm, with the greatest reduction 
occurring between 1 and 20 atm. 

The molar emission index (the moles of nitrogen oxides formed per mole of 
CO consumed) increased as the pressure was increased above 1 atm, reached a 
maximum at approximately 28 to 30 atm, and thereafter decreased slowly up to 
50 atm. 

Overall average rates of formation of nitrogen oxides were found to be 
greater than nitric oxide formation rates predicted by the conventional two- 
step Zeldovich mechanism when oxygen atoms are assumed to be in equilibrium 
with oxygen molecules, but less than those rates predicted by assuming oxygen 
atoms to be in dynamic steady state with the CO-air reactions. Circumstantial 
evidence is offered to support the proposition that the concentrations of 
nitrogen oxides in the flames are at their equilibrium values for pressures 
above about 20 atm. The case for equilibrium is not as clear at the lower 
pressures. It is suspected that reaction kinetics probably influence or pos- 
sibly control the formation of nitrogen oxides at the lower pressures. 

Comparison of the present data with data in the literature for a methane- 
air diffusion flame shows that for flames of comparable flame height (8 to 
10 mm) and pseudoequivalence ratio (0.162), the molar emission index of a CO- 
air flame is significantly greater than that of a methane-air flame. 


INTRODUCTION 

When fuels are burned with air, nitrogen oxides are typically formed from 
the interaction of nitrogen and oxygen species at the high temperatures associ- 
ated with the combustion process. Some combustion devices, such as aircraft 
turbine engines and internal combustion engines, operate over a range of pres- 
sures rather than at one single pressure. Yet, the effect of pressure on the 
formation of these nitrogen oxides NOx is not well understood. Hence, there is 
a need to determine experimentally the possible effects of combustion pressure 
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on these emissions. Such information is of interest not only because of 
the desire to know the direction and magnitude of the effects which can be 
expected, but also because such information can be used as a data base for 
the evaluation of combustion models. 

An effect of pressure on the emissions of N0 X from a laminar methane- 
air flame has already been demonstrated in reference 1 . For this flame, an 
increase in pressure above 1 atm is reported to cause first an increase and 
then a decrease in the molar emission index (the moles of N0 X produced per mole 
of fuel consumed) . The range of pressures in that investigation was from 1 to 
50 atm, with the peak emission occurring at about 9 atm. A complicating uncer- 
tainty with these results, however, is that increasingly large concentrations 
of carbon particles formed in the flame as the pressure was increased. These 
carbon particles may perturb the NO x -forming reactions from those which domi- 
nate when particulate carbon is not present. As one means of avoiding this 
complication, fuels which do not form carbon, such as carbon monoxide (CO) and 
hydrogen (H 2 ) t could be investigated over a similar range of pressures to 
determine whether the observed maximum in the molar emission index is, in fact, 
related to the presence of carbon particles. 

The present report details the results obtained from an experimental mea- 
surement of the molar emission index from a CO-air diffusion flame over the 
pressure range from 1 to 50 atm. As in reference l , the diffusion flame was 
selected for initial study of the effect of pressure on NO x formation both 
because it is representative of the type of combustion encountered in many 
practical combustion devices (aircraft turbine engines, internal combustion 
engines) and because it is much easier to study safely at high pressures than 
are premixed flames (no tendency to flashback) . This investigation was con- 
ducted in the same high-pressure facility as was used for the investigation of 
reference 1 . A qualitative explanation of the observed results is given. 


SYMBOLS 

A cross-sectional area for flow within quartz shield, cm 2 

Ab projected area of base of flame, cm 2 

Af surface area of flame, cm 2 

C-[ concentration of species i, moles/cm^ 

Cp molar heat capacity, cal/mole-K 

d thermocouple bead diameter, cm 

h convection heat-transfer coefficient, cal/cm 2 -sec-K 

1 emission index, (g of NO x as if all were N 02 )/(kg of CO consumed) 

I m molar emission index, (moles of NO x produced) / (mole of CO consumed) 
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molar emission index for N0 X at equilibrium, 

(moles of NOx produced) / (mole of CO consumed) 

equilibrium constant for oxygen dissociation reaction 

reaction rate constant for ith reaction, cm 3 /mole-sec for two-body 
reactions, cm^/mole^-sec^ for three-body reactions (number in 
subscript refers to equation number) 

radiation path length, mm 

mass flow rate of species i, g/sec 

total mass flow rate in quartz shield, g/sec 

moles of species i, moles 

molar flow rate of species i, moles/sec 

pressure, atm (1 atm = 1 01 . 3 kPa) 

universal gas constant, 82.06 cm 3 -atm/mole-K 

diffusional rate of CO into flame zone, moles/cm 3 -sec 

theoretical estimate of CO consumption rate based on global rate 
expression (see appendix D) , moles/cm 3 -sec 

reaction rate for formation of species i, moles/cm 3 -sec 

rate of formation of NO by Zeldovich mechanism, moles/cm 3 -sec 

experimental average rate of consumption of CO, moles/cm 3 -sec 

experimental average rate of formation of NO x , moles/cm 3 -sec 

theoretical estimate of CO consumption rate based on elementary 
reaction steps (see appendix D) , moles/cm 3 -sec 

temperature, K 

adiabatic flame temperature, K 
burner temperature at exit, K 
thermocouple temperature, K 
flame temperature, K 

temperature required to force agreement between I m/ e q and T m' K 
reference temperature, K 
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T w wall temperature within chamber, K 

Vrx reaction volume, cm 3 

V ST p volumetric flow rate at standard temperature and pressure, 293 K 

and 1 atm, standard cm 3 /min 

v velocity, cm/sec 

x^ mole fraction of species i 

x 0^eq Ch equilibrium oxygen-atom mole fraction for stoichiometric premixed 
combustion 

XQ, re q oxygen-atom mole fraction required by equation (4) to give the 
observed N0 X formation rate 

AHp x heat of combustion of CO, cal/g 

e emissivity 

X thermal conductivity, cal/cm-K-sec 

y viscosity, g/cm-sec 

p density, g/cm 3 

a Stef an-Boltzmann constant, 1.355 x 10“ 12 cal/cm 2 -K 4 -sec 

T residence time, msec 

cb equivalence ratio 

(}>' pseudoequivalence ratio 

Subscripts: 

act actual conditions 

eff combustion product effluent 

f flame conditions 

i ith species or ith reaction 

in flowing into flame zone 

0 ,eq based on assumption of oxygen atoms in equilibrium 

with molecular oxygen 

0,ss based on assumption of nonequilibrium steady-state balance 

for oxygen atoms in flame 
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o flowing out of flame zone 

stoich stoichiometric conditions 

Superscript: 

r reverse of reaction as written 

Chemical symbols: 

CH4 methane 

CO carbon monoxide 

CO2 carbon dioxide 

H hydrogen atom 

H2 hydrogen molecule 

H2O water 

M arbitrary third body 

N nitrogen atom 

N2 nitrogen molecule 

NO nitric oxide 

NO2 nitrogen dioxide 

N 0 X total nitrogen oxides (NO and NO2) 

0 oxygen atom 

O2 oxygen molecule 

OH hydroxyl radical 

Abbreviations: 

HC hydrocarbon 

ppm parts per million 

seem standard cm 3 /min 

i.d. inside diameter 



o.d. outside diameter 

diam diameter 

A bar over a symbol denotes time average. 


DESCRIPTION OP SYSTEM 

The high-pressure flame system used in the present study is the same as 
that described in detail in reference 1 . It consists of a diffusion flame 
burner mounted within a pressure chamber , with provisions for withdrawing com- 
bustion products and for viewing and photographing the flame. The burner is a 
3 .06-mm-diameter tube from which the fuel, CO, issues coaxially into a coflow- 
ing stream of air, which in turn is contained within a 20 . 5-mm-diameter cylin- 
drical quartz tube concentric with the burner tube. (See fig. 1.) With a 
greater than stoichiometric flow rate of air, it is possible to establish on 
this burner arrangement a confined laminar diffusion flame, referred to as a 
Burke and Schumann flame (ref. 2) . This is the typical over ventilated diffu- 
sion flame, similar in shape to a candle flame. 

To designate the overall fuel-to-air ratio, a pseudoequivalence ratio 
is defined (in keeping with ref. 1) as the actual molar ratio of fuel to air 
divided by the stoichiometric molar ratio of fuel to air. ^he term "pseudo- 
equivalence ratio" (rather than "equivalence ratio" only) is used to draw 
attention to the fact that the concept of equivalence ratio does not have the 
same significance for a diffusion flame as it does for a premixed flame. In an 
overventilated diffusion flame (i.e., 6 1 < *! ) , it is generally assumed that 

the effective average equivalence ratio in the flame zone where reaction occurs 
is, in reality, very close to unity. 

A schematic diagram of the high-pressure chamber, burner, and sample col- 
lection system is shown in figure 2. Ignition is accomplished by means of a 
resistance-heated hot-wire igniter which is retractable through a chevron pack- 
ing gland at the top of the chamber. To ignite the flame, the igniter is 
lowered to a position 1 0 to 20 mm above the burner and electrical power is 
applied. Following ignition, the igniter is retracted out of the flow of 
combustion products. The combustion products exit from the top of the quartz 
tube (or chimney) and are swept upward into a total-sample collector by the 
nitrogen gas used to pressurize the system. 

The diluted combustion products flow out of the system through a 5-ym 
particle filter and through a hot-air-heated back-pressure regulator used to 
maintain the chamber pressure at the desired level. At the exit of the back- 
pressure regulator, the sample enters a sample transfer line at near atmo- 
spheric pressure and flows to a gas analysis console. All surfaces contacted 
by the sample, from the chimney exit to the gas analysis console, are of stain- 
less steel. Color photographs of the flame are taken through safety-sight 
glass windows located on opposite sides of the high-pressure chamber. Flow 
rates of CO, air, and N 2 are measured with calibrated linear thermal mass flow- 
meters having an accuracy of 1 percent of full scale. Meter ranges are 0 to 
200 seem for CO, 0 to 3000 seem for air, and 0 to 50 000 seem for N 2 . 
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The high-pressure flame system is operated by setting the desired chamber 
pressure with the back-pressure regulator and pressurizing with N 2 - After 
reaching pressure, flow rates for CO, air, and N 2 are adjusted to their target 
values: 73 seem or 146 seem for CO, the desired flow rate for air, and about 
6000 seem for N 2 . The flame is then ignited. Two different cylinders of CO 
were used in the experiments: one contained 99.3 ppm H 2 with 327 ppm total 

hydrocarbon (HC) contamination; the second contained less than 10 ppm H 2 with 
less than 4 ppm HC contamination. The trace amounts of H 2 were intentional; 
their purpose was to serve as a chain carrier in initiating and sustaining the 
combustion of the CO. No difference was observed in the data obtained from 
either of these two CO cylinders. The air was zero grade, 99.995 percent pure 
with an HC concentration of less than 1 ppm; the N 2 was also 99.995 percent 
pure. To prevent both the appearance of a reddish-yellow luminosity near the 
tip of the flame and a reddish-brown deposit on the burner exit rim (apparently 
caused by the presence of trace amounts of iron carbonyl in the CO) , the CO 
from the supply cylinder was first passed through a dry-ice trap and then 
through an adsorbent bed of Molecular Sieve 5A. 


EXPERIMENTAL MEASUREMENTS 

Experimental data consist of the chemical composition of the combustion 
products, color photographs of the flame, and flame temperatures. 


Chemical Analyses 

Chemical analyses for CO 2 , CO, 02, and N0 X were performed on the total gas 
mixture exiting the high-pressure chamber. Analyses for CO and CO 2 were per- 
formed on separate (commercial) nondispersive infrared analyzers, and the anal- 
ysis for O 2 was performed on a (commercial) paramagnetic oxygen analyzer. For 
the measurement of NO x , a slightly modified, commercially available chemilumi- 
nescent NO— NO x analyzer was used. 

Calibration gases for the infrared and paramagnetic analyzers were gas 
mixtures of CO, CO 2 , and O 2 in N 2 , with a stated analytical accuracy of better 
than ±2 percent. The calibration gases for the NO x analyzer were gas mixtures 
of NO in N 2 , with a stated accuracy of ±1 percent. During an experimental run, 
readings of the calibration and zero gases were frequently interspersed between 
readings of the sample gas to eliminate possible errors caused by instrument 
drift. 


Photographs 

Color photographs of the flame were taken with a 35-mm single-lens-reflex 
camera fitted with a 200-mm/f2.8 reflective lens and with color negative film. 
Since the intensity of the bluish color of the flame was observed to change 
very little with pressure, a constant camera exposure of 1 /I 5 sec was used at 
all flame conditions. 


7 


1 



Flame Temperature 


Because of the difficulty of gaining physical access to the flame and 
the steep thermal gradients in the flame, probing with thermocouples presents 
severe experimental difficulties and was not attempted. Optical techniques 
relying on incandescent carbon particles in the flame, like the one employed 
with the methane-air diffusion flame in reference 1, could not be used because 
of the absence of such particles. Sodium line reversal was ruled out because 
introduction of sodium into the flame presents certain experimental difficul- 
ties: the presence of sodium in the system could perturb the N0 X formation 

rates and the deposition of sodium throughout the chamber and sample lines was 
to be avoided. These practical considerations led to the decision to settle 
for estimating the flame temperature at all pressures from a knowledge of the 
flame temperatures measured at 1 atm. To this end, the burner was removed from 
the high-pressure chamber and was set up on a laboratory bench at ambient pres- 
sure. Flames were established on the burner at the same flow rates of fuel and 
air as in the high-pressure experiments, and were probed with a chromel-alumel 
thermocouple with a wire diameter of 0.076 mm and a bead diameter of 0.269 mm. 
No attempts were made to obtain temperature profiles throughout the flames, and 
only peak temperatures were recorded. 


RESULTS 

Experimental measurements were made at two flow rates of CO over the pres- 
sure range from 1 to 50 atm, 73 seem and 146 seem. The limits of stable burn- 
ing and the shapes of the flames changed markedly with increasing pressure. In 
this section, the results of the measurements of stability and extinction lim- 
its, flame shape, NO x , and flame temperature as functions of pressure are pre- 
sented. The N0 X data are treated in the section "Discussion." 


Stability and Extinction Limits 

A rough determination of the stability and extinction limits was made with 
the burner prior to selecting the final experimental CO and air flow rate set- 
tings. These stability and extinction limits were determined at l atm, 20 atm, 
and 50 atm; the results are shown in figure 3. Much of the scatter in the data 
was caused by an air line leak which was discovered after the data had been 
taken; the data were then corrected to account for this leak. Hence, the indi- 
cated limits are only approximate. However, the important feature to observe 
in figure 3 is that the fuel-air flow rate regimes in which combustion occurs 
decrease considerably as pressure increases. This phenomenon is similar to 
that reported in reference 1 for a methane-air diffusion flame on the same 
burner system. 

On the basis of the data in figure 3, two fuel flow rates were selected 
for study: 73-sccm and 1 46-sccm CO. Corresponding to these fuel flow rates, 

flow rates of air were selected at 1 730 seem and 21 40 seem, respectively, giv- 
ing pseudoequivalence ratios <b f of 0.100 and 0.162, respectively. At 73-sccm 
CO, the air flow rate of 1730 seem is well within the stable region while at 
1 46-sccm CO, the air flow rate of 21 40 seem is near the extinction boundary at 
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the higher pressures. A pseudoequivalence ratio 4>* of 0.100 could not be 
easily achieved at the higher fuel flow ratd’ of 1 46-sccm CO because an air 
flow rate of 3476 seem would be required, a rate above the limit of the air 
flow meter. 


Flame Shape 

Marked changes in shape and structure of the CO-air diffusion flames occur 
as the pressure is increased from 1 to 50 atm. These changes can be easily 
seen in the color photographs in figure 4. The bottom edge of each photograph 
is coincident with the top edge of the burner, and the lateral position of the 
burner is shown below each photograph. The bluish color of the flames is prob- 
ably caused by O 2 Schumann-Runge radiation and CO + O continuum radiation. 
(See ref. 3.) 

From 1 atm to about 5 atm, the flames bow outward from the inner fuel core 
into the surrounding air annulus, with the result that the diameter of the 
flames is larger than the diameter of the burner exit. At about 10 atm, the 
diameter of the flames is reduced, and the shapes of the flames closely approx- 
imate cones. At pressures above 10 atm, there is an increasing concavity to 
the flames, and their diameters decrease noticeably. A similar behavior was 
observed for the methane-air diffusion flame in reference 1 , wherein it was 
conjectured that the change in shape might result from the presence of carbon 
particles in the flame. In view of the present results, this conjecture is 
obviously incorrect, and some other explanation must be found. (It should be 
noted that, according to classical Burke and Schumann flame theory, no change 
in the size or shape of a confined laminar diffusion flame should occur as the 
pressure is changed.) Comparison of the shapes of the present flames with 
those of the methane-air flames discussed in reference 1 suggests that the only 
effect which can be attributed to the presence of carbon particles is a slight 
bulge occurring in the upper third of the methane-air flames somewhat below 
their tips. 

The heights of the CO-air flames decrease with pressure from 1 atm to 
about 1 0 atm and thereafter remain relatively constant up to 50 atm. These 
flame heights are plotted in figure 5. The decrease in the height of the flame 

1 1 

at a 1 46-sccm flow rate is from about 9— mm to about 7— mm, or roughly 20 per- 

2 2 

cent; the decrease for the flame at a 73-sccm flow rate is from about 6 mm to 
about 3 mm, or roughly 50 percent. 


Oxides of Nitrogen 

Results of the experimental measurements on the total oxides of nitrogen 
are shown in figure 6 as I m , molar emission index of NO x , as a function of 
pressure. Experimentally, this index is determined as the ratio of the mea- 
sured concentrations of NO x to CO 2 in the effluent from the high-pressure cham- 
ber. The molar emission index is analogous to the emission index customarily 
used to report data on gas turbine engines; the latter is defined as the number 
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of grams of NOj^ produced, as if all nitrogen oxides were present as NO 2 , per 
kilogram of fuel burned. According to these definitions, the emission index 
for a CO flame is 1 045 times the molar emission index. 

To obtain the data points in figure 6 , the CO flow rates were set at their 
target values, three or more runs were made at different air flow rates, and 
concentrations of NOjjr CO, CO 2 , and O 2 were measured. Then values of I m cal- 
culated from these data were plotted against pseudoequivalence ratio <b 1 and 
were read back at the desired values of (f ) 1 (appendix A gives the method of 
calculating (t* 1 ). This procedure effectively eliminated the need for setting 
the air flow rates at precise values, a difficult task at times. It also 
avoids a reliance on the accuracy of the air flow meter. On the whole, over 
the usual ranges of (J) ' employed (roughly 0.08 to 0.14 for 73-sccm CO and 0.14 
to 0.18 for 1 46-sccm CO), there was no pronounced or consistent variation of 
molar emission index with 4) 1 which could be attributed to other than experi- 
mental error. Above about 30 atm, however, there did appear to be a general 
decrease of I m with <t>', although this decrease was not large. For purposes 
of comparison, the molar emission index for the methane-air flame studied in 
reference 1 has been included in figure 6 . 

Figure 6 shows that for both CO flames, the molar emission index starts 
at a relatively low value at low pressure and increases as the pressure is 
increased. The molar emission index for the flame with the lower flow rate 
reaches a maximum of about 2.5 x 1 0“ 3 , while the I m for the flame with higher 
flow rate reaches a maximum of about 4.6 x 10“^, both maxima occurring at 
approximately 28 to 30 atm. At still higher pressures, the I m for both 
flames levels off and drops slightly. These results are in contrast with those 
for the CH 4 flames, which exhibit a much smaller increase with pressure, attain 
an earlier maximum (at about 9 atm), and then decrease to near the low-pressure 
value again at 50 atm. This comparison between the GH 4 flame, with its low 
flow rate of 42 seem, and the CO flames with a higher flow rate may at first not 
seem proper; however, it is noted that in terms of flame height (8 to 1 0 mm), 
pseudoequivalence ratio (0.162), and air flow rate (2450 seem), the CH 4 flame 
is quite similar to the 1 46-sccm CO flame. Yet the I m for this CO flame is 

1 

roughly 4— times the I m for the CH 4 flame. One obvious, at least partial, 

explanation is that the CO flames probably attain higher flame temperatures by 
virtue of their higher adiabatic flame temperatures - 2384 K for CO as compared 
with 2230 K for CH 4 J Furthermore, because radiation from the CH 4 flame is 
significant, its actual flame temperature is likely to be even lower than the 
1 54 K difference in adiabatic flame temperatures would indicate. At the lower 
pressures, where carbon radiation from the CH 4 flames is not yet substantial, 
the I m values for both flames are in closer agreement. The effects of tem- 
perature on I m are covered in more detail in the section "Discussion.” 


^ At 1 atm. Although there is a slight increase in adiabatic flame temper- 
ature with increase in pressure for both fuels, the increase is not dissimilar. 
The values of adiabatic flame temperature were calculated using the chemical 
equilibrium program of Svehla and McBride (ref. 4) . 



Only two studies could be found in the literature where measurements of 
N0 X in CO-air flames at elevated pressures are reported. These data are com- 
pared with the present data in figure 7 on a log-log plot. The data points A, 
B, and C are from Haber and Coates (ref. 5) for a diffusion flame burner simi- 
lar in design to the one used in the present study f but with a 1 . 5-mm-diameter 
fuel port. Most of the Haber and Coates data are not for air but for support- 
ing atmospheres consisting of O 2 /N 2 mixtures of varying composition. Points A, 
B, and C are the only data points at which direct comparisons with the present 
results are possible. Even though these data are for different flow rates of 
CO, and were obtained with a burner of smaller diameter, they are within rea- 
sonable proximity of the present data. 

The data of Newitt and Lamont (ref. 6), shown as point D, are for an 
opposed-jet diffusion flame and represent their only data for air. Although 
they did not state the flow rates of their reactants, it is reasonable to 
assume that because of the opposed-jet flow, there must have been a turbulent, 
fairly well-mixed reaction zone. This type of reaction system might be 
expected to approach more closely that of premixed combustion rather than 
diffusional combustion, or, at the very least, to be somewhere between the two. 
However, there are other problems with the data of Newitt and Lamont that make 
an interpretation of their data difficult. In particular, large concentrations 
of unburned CO remained in the product gases even at relatively low values of 
<b', demonstrating that combustion was far from complete. 


Flame Temperature 

As mentioned previously, temperature data were obtained for flames at both 
flow rates at 1 atm only. With these data and estimates of heat losses from 
the flames, a reasonable trend of flame temperature with increasing pressure 
was calculated. The raw thermocouple data were first corrected for radiation 
losses by the methods detailed in appendix B. Flame temperature as a function 
of pressure was then estimated by the procedure described in appendix C. No 
attempt was made to correct the thermocouple data for the heat conduction 
losses that can arise from the steep temperature gradients in the flame. (In 
severe cases, such losses could possibly produce indicated readings as much as 
1 00 K too low, but these are extremely difficult to estimate accurately.) The 
resulting temperature curves are shown in figure 8 along with a plot of the 
adiabatic flame temperature for a CO-air flame at d> = 1.0. 2 These temperature 
curves are believed to be fairly reasonable estimates of the true flame temper- 
atures. For the i 46-sccm CO flame, temperature increases approximately 370 K 
between 1 atm and 50 atm; for the 73-sccm CO flame, there is a corresponding 
increase of approximately 480 K. These increases compare with an increase of 
1 40 K in adiabatic flame temperature over the same pressure range. 


2 It may bear repeating here that even though the overall, or pseudoequiva- 
lence ratio 6' may be considerably smaller than unity, the effective average 
equivalence ratio within the reaction zone must be quite close to unity. The 
adiabatic flame temperature was calculated from the equilibrium computer pro- 
gram of Svehla and McBride (ref. 4) . 


n 



DISCUSSION 


In this section, the average rate of formation of N0 X is estimated from 
the molar emission index and the average rate of consumption of CO. The ques- 
tion of the combustion reactions approaching equilibrium is treated, and the 
average rates of formation of NOx are discussed and compared with various theo- 
retical predictions. Finally, the shapes of the molar emission index curves 
are discussed from the point of view of their being influenced by approach to 
equilibrium, particularly at the higher pressures. 

Average Rate of Formation of N0 X 

An estimate of the average rate of formation of N0 X within the flame was 
obtained as follows. At steady-state conditions and after complete combustion 
of the fuel. 


I 


m 


n N0 x 


n co 


. • 

where nro is the molar flow rate of CO into the flame and njjo v is the 

molar flow rate of N0 X out of the flame. After rearrangement of this expres- 
sion and division by a reaction volume V Rx , there results 


• • 

n N0 x n C0 

V Rx " "* V Rx 

If the reaction volume in which NC>x is formed is taken to be approximately 
equal to (though not necessarily coincident with) the reaction volume in which 
CO is consumed, this expression can be rewritten as 

— E — E 

r N0 x * *in r CO H > 


where r^Q is the experimental average rate of formation of N0 X and r^ is 

the experimental average rate of consumption of CO. An estimate of the magni- 

— E * 

tude of r^o was obtained by dividing the known molar flow rate of CO n co by 

the volume of the flame zone, experimentally obtainable from the color photo- 
graphs in figure 4. (These reaction volumes should be a fairly accurate mea- 
sure of the volumes in which fuel is actually consumed since the blue zones in 
figure 4 are largely the result of chemiluminescence from the CO + 0 reac- 
tion, and since combustion was observed to be virtually complete within the 

— E 

flame.) Calculation of r N Q x from equation (1) implies that the reaction 
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volume for N0 X formation is roughly equal to that for CO combustion. Although 
this is probably not strictly true, it is certainly reasonable to assume that 
these two volumes are roughly proportional. Since these reaction volumes prob- 
ably do not differ by much more than a factor of 2 either way, in the absence 
of information to the contrary, they have been taken, for convenience, to be 
equal. 


The photographs in figure 4 were used to determine the magnitudes of the 
CO reaction volumes by tracing the thin blue reaction zones at the sur- 

faces of the flames and by analytically revolving the areas so defined about 
the center line of the flames. The resulting reaction volumes are shown in 

— E 

figures 9(a) and 9(b). Values for r G o based on these reaction volumes 
are plotted as a function of pressure in figure 10(a). Also plotted in fig- 
ure 10(a) are two different predictions of CO combustion rates r G and r T 
for a global rate expression (ref. 7) and a set of elementary reaction steps 
for CO combustion, respectively. Both are based on the assumption of an effec- 
tive average equivalence ratio of unity in the flame zone and on initial com- 
position. Details of the methods used for evaluating r G and r T are given 
in appendix D. 

Because r q and r«p are based on initial composition and an average c|) 
of unity, they very likely represent the approximate maximum rates of combus- 
tion of CO in air for the present flame conditions. Departure of the experi- 

— E 

mental rates Tqq from these predicted rates could mean that the assumption of 
an average (f> of unity is not a good one, that the choice of initial composi- 
tion for the reactants is not appropriate, that the use of Tf as an average 
temperature throughout the flame zone is not proper, that diffusion is really 
controlling the rate of reaction, or some combination of these. To treat the 
reaction zone in a diffusion flame in such simplified terms as these obviously 
oversimplifies the true situation. However, this simplified picture allows the 
inference from figure 10(a) (by virtue of the fact that the predicted rates 
are, for the most part, greater than the experimental rates) that the rates of 
combustion of CO in the present flames are largely diffusion controlled, cer- 
tainly at pressures above 2 atm or so. 

An additional reason to suspect diffusional control of the present CO-air 

— E 

flames, even at all pressures, is the proximity of the observed rates r G0 
for flames at both flow rates. Because the flames have significantly dif- 
ferent flame temperatures they also would be expected to have much different 
reaction rates if reaction kinetics were actually controlling. However, since 
these rates do not differ significantly, and since rates of diffusion are rela- 
tively insensitive to temperature, it is reasonable to infer that the combus- 
tion rates in both flames are largely diffusion controlled. A stronger argu- 
ment for diffusional control can be obtained by estimating the diffusional 
rates of CO into the flame zone r D . This estimate was carried out by the 
method described in appendix E, and the resulting values of r D are shown in 

-E 

figure 10(b). (Because these r D virtually overlay the r^Q in figure 10(a) 
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they have been plotted on a separate figure for ease in visualization.) The 

— E 

almost exact agreement between r co and the estimated rates r D seems con- 
vincing evidence that the combustion reactions are, in fact, diffusion con- 
trolled. The important implication of the flames being diffusion controlled 
is that the combustion reactions approach equilibrium. 


The average rates of formation of nitrogen oxides in the flame were esti- 

— E 

mated using equation (1) and figures 6 and 10(a) for I m and tqq, respec- 

— E 

tively. The resulting rates r N Q x are plotted in figure 11. Also plotted in 

figure 11 are two sets of theoretical rate predictions, each based on a differ- 
ent set of assumptions, but both employing the conventional two-step Zeldovich 
mechanism for NO formation: 


O + N 2 * NO + N 


( 2 ) 


N + 0 2 * NO + O (3) 

By standard treatments (see, for instance, ref. 1), the rate of formation of NO 
by the Zeldovich mechanism is given by 


r NO = 2k 2 c O C N 2 


(4) 


where k 2 , the forward reaction rate constant of reaction (2), is (refs. 1, 8, 
and 9) 


k 2 = 1.36 x TO 14 e -37947 / 1. 

One set of the theoretical rate predictions in figure 11 is based on the 
assumption that oxygen atoms are in equilibrium with molecular oxygen, and the 
other is based on the assumption that there exists superequilibrium steady- 
state concentrations of oxygen in the flames. To obtain the theoretical rate 
predictions for the situation of oxygen atoms in equilibrium with molecular 
oxygen, the equilibrium constant (refs. 1 and 9) 


Kq = 


(Co ) 2 


3.56 x 10^ — e -59386/T 
T 


is substituted into equation (4) to eliminate the unknown oxygen-atom concen- 
tration Cq. In terms of mole fractions, this rate expression then becomes 


r N0 = 6.90 x 1 o 1 3 


p3/2 

T 2 


e -67640/T 



x N 2 


(5) 



This expression was evaluated at the maximum possible values Xq 2 = 0.21 and 

Xfl 2 “ 0-79, at T = Tf, and is labeled ( r N0^0,eq* Rates so obtained are the 

maximum rates of formation of NOx attainable when oxygen atoms are in equilib- 
rium with molecular oxygen since the effective mole fraction for oxygen in the 
zone where nitric oxide is formed must be substantially less than the 0.21 
used, and the effective mole fraction for nitrogen is probably somewhat less 
than the 0.79 used. For instance, in reference 1 it was estimated that in a 
methane-air diffusion flame, a reasonable value for the effective oxygen mole 
fraction is on the order of 0.04. If this were also the case for the present 
CO-air diffusion flames, the predicted rates (^No)o,eq fi 9 ure would be 
about 0.4 times those indicated - lower yet, and farther away from agreement 

— E 

with the experimentally observed average rates r^Ox* Even though these exper- 
imentally observed rates are rough estimates, it seems unlikely that they could 
be in error by a factor as large as 100 or more, which they would have to be at 
1 atm if the predicted rates ( r NO)o,eq are correct. 

One of the more reasonable explanations for the higher experimental rates 
is that the oxygen atoms in the flame zone are not in equilibrium with the 
molecular oxygen but are present in concentrations greatly in excess of equi- 
librium values. Such an explanation is consistent with the observation that 
the experimental rates approach the predicted rates ( r NO^O,eq at tlle higher 
pressures where superequilibrium oxygen-atom concentrations are less likely. 
Calculations of the magnitude required for these superequi librium oxygen-atom 

— E 

concentrations for the predicted rates to agree with r N Q x can be made by sub- 
— E 

stituting r NO x for r NO equation (4), taking t = Tf and x N2 = 0.79, and 

solving for the oxygen-atom concentration. The resulting concentrations are 

stoich 

plotted in figure 1 2 as the ratio x o,req/ x O,eq ' where the denominator is the 
equilibrium oxygen-atom mole fraction for stoichiometric premixed combustion at 
Tf. Extremely high values of the required oxygen-atom ratio occur at the lower 
pressures, but they decrease continually as the pressure is increased. 

To obtain a more realistic estimate of the oxygen-atom concentrations, an 
oxygen-atom balance was made. This balance is based on the assumption that a 
steady-state condition exists in the flame wherein the rate of generation of 
oxygen atoms is exactly counterbalanced by the rate of removal of oxygen atoms. 
The details of this approach are given in appendix F. By this method oxygen- 
atom concentrations have been calculated and converted to the ratio 

s to i ch 

x O,ss/ x O,eq and are plotted in figure 12. For the 1 46-sccm CO flame, this 
concentration of oxygen atoms is considerably greater than that required to 
produce the experimentally observed rates, and for the 73-sccm CO flame, this 
concentration is greater than that required for all pressures above about 
3 atm. For pressures below about 3 atm, the estimated oxygen-atom concentra- 
tion is less than that required, being as much as an order of magnitude less at 
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1 atm. This lower value of xo fSS for the 73-sccm CO flame below 3 atm implies 
that there is an insufficient concentration of oxygen atoms generated by the 

combustion process below 3 atm to account for the observed rate r^Ov b Y the 

Zeldovich mechanism. However, because of the size of the potential errors 

—E 

involved in calculating rjjo v f i- n selecting an appropriate value for x N , in 

making the steady-state analysis for xo fSS , as well as in estimating Tf, there 
is insufficient cause, on these grounds alone, to invalidate the Zeldovich mech- 
anism as being the basic mechanism for the production of nitric oxide. 

Rates of formation of nitric oxide based on these steady-state oxygen-atom 
concentrations (cno)o,ss have been calculated using equation (4) and are com- 

— E 

pared with the experimentally observed rates r^Ov figure 11. For this cal- 

culation, the nitrogen concentration was taken to be - 0.556 (i.e. f that 

nitrogen concentration existing in a stoichiometric mixture of CO and air) . 

These rates, ( r NO)o,ss' are not onl Y considerably greater than the rates pre- 
dicted on the basis of oxygen-atom equilibrium but are also greater than the 
experimentally observed rates except for the 73-sccm CO flame below about 3 atm. 
Clearly, the reason why (r^o^O S s * s considerably greater than UNo)o,eq * s 
that the oxygen-atom concentrations based on the steady-state analysis are con- 
siderably greater than those based on the equilibrium assumption. (Examination 
of the analysis in appendix F shows that this is the case because the reaction 
step with the fastest rate of production of oxygen atoms is the combustion 
reaction 


CO + 0 2 + C0 2 + 0 
and not the dissociation reaction 
0 2 + M + 20 + M 

which it would have to be for atomic oxygen to be in equilibrium with molecular 
oxygen . ) 

From figure 11, it is apparent that if the Zeldovich mechanism is the cor- 
rect mechanism for nitric oxide formation in CO-air diffusion flames, and if 

“E 

the experimentally determined rates r^o are approximately correct, then the 

true average oxygen-atom concentrations in the flames are not as large as those 
calculated on the basis of the steady-state analysis, that is, on the basis of 
purely kinetic considerations. This is certainly the case at pressures above 

-E 

about 15 atm or so, where (ruo)o ss is sufficiently greater than r^v so 
that it is unlikely that the difference can be attributed to experimental 
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error. Below 15 atm, and particularly below about 3 or 4 atm, the case is not 
so clear. Precisely the same conclusion can be reached by inspection of fig- 
ure 12, where x o, re<j/ x 0^eq Ch is less than x 0 S s/ x 0^eq Ch at a11 except the 
lower pressures. This result leads one to consider the possibility that per- 
haps reaction kinetics may not be controlling the rate of formation of nitric 
oxide, but that diffusion might be playing a role. If diffusion were, in fact, 
controlling, then it necessarily follows that the NOx concentrations in the 
flames should approach their equilibrium values. This eventuality is consid- 
ered in the following section. 


Molar Emission Index 

Equilibrium values of the molar emission index calculated for stoichiomet- 
ric combustion and T = Tf are plotted in figure 1 3 as I m ,eq along with the 
curves of l m from figure 6. The similarity in shape between these curves is 
obvious. However, I m ,eq rises faster and peaks sooner than does I m . Above 
about 28 atm, the curves of I mf eq parallel those of I m fairly well but are 
roughly 20 to 25 percent lower. Yet, the general similarity in shape suggests 
the possibility that the NO x concentrations in the flames might be approaching 
their equilibrium values, particularly at the higher pressures where the curves 
of I m and I mf eq are reasonably parallel, and where the required oxygen-atom 

ratios x o, req/ x 0^eq C ^ *- n figure 12 are decreasing to lower, not so excessive 
values. At pressures above 28 atm, separate calculations indicate that an 
increase in Tf of only about 50 K for both flames would be required to bring 
I m ,eq into exact agreement with I m , a variation certainly well within the 
error of estimating Tf. Circumstantial evidence of a similar nature (i.e., 
the agreement of Im,eq and x m) was offered in reference 1 for a CH 4 -air dif- 
fusion flame in support of the proposition that the NO x concentrations in that 
flame had approached equilibrium levels at pressures above about 20 atm or so. 
Unlike the present CO-air flames, however, it was clear for the CH 4 -air flame 
that equilibrium levels of NOx had not been approached at the lower pressures, 
and that kinetics were still controlling. For the present CO-air flames 
the case is not so clear. However, because of the very large values of 

x 9, reg/ x 0^eq Ch at the lower pressures (see fig. 12), it is suspected that 
kinetics are controlling even though lm,eq reasonably close to I m . 

Because the experimentally estimated temperatures Tf are approximate in 

nature and because equilibrium levels of nitric oxide are quite sensitive to 

temperature, additional calculations were made to determine those temperatures 

which would have had to exist in the flames in order for I m ,eq to coincide 

identically with I m over the entire range of pressures. These temperatures 

are plotted in figure 1 4 as T T . Comparison of these temperatures with Tir 

A m J - 

shows the average deviation between them to be about 40 K, a value well within 
the experimental error of the estimated temperatures Tf. Such close agreement 
between T T and Tf is, nevertheless, circumstantial and therefore insuffi- 

cient to prove that the NO x concentrations approach their equilibrium values at 
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all pressures* Further, because of the uncertain nature of Tf and the great 
sensitivity of I m , eq to T f' case for or against equilibrium N0 X cannot 
be resolved simply by consideration of the molar emission index. At the lower 

pressures, the very large values of xq^ re g/xQ^ eq° h tend to argue against 
equilibrium there. Yet, the possibility that NOx concentrations approach equi- 
librium within the flames at the higher pressures above about 20 atm seems dis- 
tinctly real. 

To rule out the possibility that the shape of the I m curves might be 
caused by changes in residence time, residence times were calculated from the 
expression 


T 


Vrx 

V4>=r1 


where Vp* is the reaction volume (plotted in figs. 9(a) and 9(b)) and 
is the volumetric flow rate of CO and air at an average equivalence ratio of 
unity and at T = Tf. Resulting values of x are shown in figure 15, where x 
first decreases, then passes through a minimum between 5 and 7 atm, and there- 
after continues to increase with increasing pressure. Since this trend is con- 
trary to the one for I m , it is obvious that the shape of the I m curves is 
not determined by residence-time considerations; in particular, the decrease in 
I m observed at the higher pressures is not caused by a corresponding decrease 
in X. 


To help resolve this question of the NOx-forming reactions approaching 
equilibrium, a mathematical model of a CO-air diffusion flame would be 
extremely useful. Tf such a model employed reasonable kinetic mechanisms 
describing combustion, and if the model adequately predicted the observed 
changes in flame structure with pressure, then there would be some hope that 
the basic details of the NO x formation processes could also be adequately 
modeled, and, hence, better understood. Recently, a variety of analytical 
approaches and computer codes capable of mathematically modeling both the com- 
plicated fluid dynamic and chemical aspects of diffusion flames have become 
available (refs. TO to 19), and these could possibly prove useful. 

As a final comment, it was speculated in reference 1 that among the vari- 
ous factors which could possibly account for the presence of the maximum in I m 
for a CH 4 ~air diffusion flame is the large amount of carbon in the flame. The 
present results clearly show that the presence of carbon in the flame is not a 
requirement for such a maximum. Although such carbon could still conceivably 
affect the location of the maximum, it is obviously not necessary for such a 
maximum to exist. Based on the present data and analysis, the preferred 
explanation for the shape of the X m curves as a function of pressure for both 
the CO-air and Cl^-air diffusion flames is that equilibrium levels of NO x are 
approached in the flames at the higher pressures. 



CONCLUDING REMARKS 


A study has been made of nitric oxide formation in a laminar CO-air diffu- 
sion flame over a range of pressures from 1 to 50 atm. The carbon monoxide 
(fuel) issued from a 3. 06-mm-diameter port coaxially into a coflowing stream of 
air confined within a 20. 5-mm-diameter chimney. Experimental measurements of 
nitric oxide concentrations were made at two carbon monoxide (fuel) flow rates: 
73 standard cm 3 /min (seem) and 1 46 seem. The flow rates were held constant 
over the entire range of pressures from 1 to 50 atm. The stability and extinc- 
tion limits, flame shape, flame height, NO^, and flame temperature were all 
affected by changes in pressure. Based on these experimental data and their 
analysis, the following observations and conclusions are offered: 

1 . The shape of CO-air diffusion flames changes from wide and convex at 
1 atm to slender and concave at 50 atm. The most noticeable change in shape 
occurs below about l 0 atm. 

2. The height of the flames decreases with pressure from 1 atm to about 

1 0 atm or so, and thereafter remains relatively constant up to 50 atm. For the 

1 1 

flame with 1 46-sccm flow rate, the decrease is from about 9— mm to about 7— mm; 

2 2 

for the 73-sccm flame, the decrease is from about 6 mm to a little over 3 mm. 

3. The regimes of stable burning decrease as the pressure is increased, 
with the greatest reduction occurring between 1 atm and 20 atm. 

4. The concentrations of nitrogen oxides produced by the flames were mea- 
sured and reported as molar emission index (the moles of nitrogen oxides formed 
per mole of carbon monoxide consumed) . For both flames, the molar emission 
index increases as the pressure is increased above 1 atm, reaches a maximum in 
the neighborhood of 28 to 30 atm, and thereafter decreases slowly up to 50 atm. 
Hence, at all pressures above 1 atm and up to at least 50 atm, the molar emis- 
sion index is greater than that at 1 atm. 

5. Comparison of the present data with data in the literature for a 
methane-air diffusion flame shows that for flames of comparable flame height 
(8 to 1 0 mm) and pseudoequivalence ratio (0.162), the molar emission index of 

1 

a CO-air flame can be as much as 4— times that of a methane-air flame. 

2 

6. Those few data in the literature on nitric oxide concentrations gen- 
erated by coaxial CO-air diffusion flames are in reasonable agreement with the 
present data. 

7. Overall average reaction rates of carbon monoxide were determined for 
the diffusion flames. Comparison of these rates with kinetic predictions and 
diffusion rates suggests that the carbon monoxide combustion reactions are dif- 
fusion controlled. The implication of this is that these combustion reactions 
approach equilibrium. 


1 9 


8. Overall average rates of formation of nitrogen oxides were estimated 
from the data. These rates are greater than nitric oxide formation rates pre- 
dicted by the conventional two-step Zeldovich mechanism when oxygen atoms are 
assumed to be in equilibrium with oxygen molecules, but are, in general, less 
than the predicted rates when oxygen atoms are assumed to be in dynamic steady 
state with the CO-air reactions. 

9. Circumstantial evidence is offered in support of the proposition that 
the concentrations of nitrogen oxides in the CO-air flames are near their equi- 
librium values for pressures above about 20 atm. The case for equilibrium is 
not clear at the lower pressures, but it is suspected that kinetics probably 
influence or even control the formation of nitrogen oxides at the lower 
pressures. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
August 1 0, 1 979 
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APPENDIX A 


METHOD OF CALCULATING PSEUDOEQUIVALENCE RATIO 

The pseudoequivalence ratio d> * was obtained from the composition of the 
effluent from the high-pressure chamber. It was calculated from its definition 


_ ^ n C 0 / n 0 2 ) act _ (X 00 /X 02 ) act 

(n C o/ n 0 2 ) stoich < X C0/ X 0 2 > stoich 

where 


x co\ 

X °2/stoich 


2 


(AT) 


The ratio of CO to O 2 in the numerator was obtained from an analysis of the 
effluents as the ratio of the unburned CO plus the burned CO (CO 2 ) to the sum 
of the excess oxygen remaining after combustion and the amount of oxygen con- 
sumed during the combust ion, as indicated in the following equation: 


x co' 

V X °2/act 


X C0 + X C0 2 


v x o 2 + r x co : 


eff 


The air flow rate was calculated from 1 
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APPENDIX B 


THERMOCOUPLE CORRECTIONS FOR RADIATION LOSSES 

To correct the raw thermocouple readings for radiation losses, the follow- 
ing formula of Kaskan from reference 20 was used: 


T f - T c = 


i . 25ecrd 3 / 4 /V N 1 /4 


/y V 7 

{— (T c 4 - T w 4 ) 

\Pv/ 


The value of e was taken to be 0.25 (refs. 21 and 22), d was obtained from 
micrometer measurement of the thermocouple bead, the product pv was obtained 
from the expression pv = m«p/A, and X and y were evaluated at Tf for a 
CO 2 and N 2 mixture resulting from the stoichiometric combustion of CO and air. 
Mixture properties for y and X were calculated using Wilke’s equation 
(ref. 23) and Brokaw's equation (ref. 23), respectively. The pure-gas proper- 
ties for y and X were obtained from reference 24. The quantity T w 4 was 
considered negligible compared with T c 4 , and was therefore dropped from the 
equation. 
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APPENDIX C 


ESTIMATION OF FLAME TEMPERATURE AS A FUNCTION OF PRESSURE 

An approximate estimate of flame temperature was obtained by mathemati- 
cally accounting for those sources of heat loss from the flame which reduce its 
temperature to some value below the adiabatic flame temperature, as follows. 

When there is no heat loss from the flame. 

Heat Sensible Sensible 

generated + heat in = heat out 

where the heat generated is from the heat of combustion and the sensible heat 
refers to the heat content of the gases flowing into the flame zone (CO and 
air) and out of the flame zone (CO 2 and N 2 ) - Stated mathematically, this 
becomes 


AHRxnicO + l Cp, in lH(T in - T ref ) = l Cp f0 A 0 (T ad - T ref ) 


(Cl ) 


When there is heat loss from the flame. 

Heat Sensible Sensible Heat radiated Heat radiated Heat convected 

generated + heat in = heat out + to chamber + to burner + to burner 

or 

^Rx m CO + I c p, in n in^ T in ~ T ref) = I c p, o n o( T f “ T ref) + Af£fCf(Tf 4 - T w 4 ) 

+ A b £fCT (Tf 4 - T b 4 ) + hA b (Tf - T b ) (C2) 

Equating equations (Cl) and (C2), collecting terms, neglecting T w 4 and T b 4 
relative to Tf 4 , and approximating T b by T im gives 

I c p,on 0 ( T ad " T f) = (A f + A b ) e f aT f 4 + hA b (T f - T in ) (C3) 

The term £ c p, o n o was evaluated for a CO 2 and N 2 mixture resulting from the 
stoichiometric combustion of CO and air and from the known flow rate for CO. 

The quantities Af and A^ were obtained from measurements of the external 
boundaries of the luminous zones in the color photographs in figure 4 and are 
shown plotted in figures 16 and 17. The emissivity Gf was obtained using the 
method of reference 25. For this purpose, the luminous zones in figure 4 were 
taken to consist of a CO 2 and N 2 mixture resulting from the stoichiometric com- 
bustion of CO and air, the radiation path length L was obtained from measure- 
ments of the thickness of these luminous zones in the upper part of the flames 
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(L is shown plotted in fig. 18), and corrections for pressure broadening 
of the CO 2 bands were made by employing correction factors obtained from 
reference 25. 

At 1 atm, where Tf is known from direct thermocouple measurement, 
equation (C3) was solved for the heat-transfer coefficient h, giving 
0.362 cal/cm 2 -min-K for the i 46-sccm CO flame, and 0.331 cal/cm 2 -min-K for 
the 73-sccm CO flame. Because the heat-transfer coefficient is a function 
of the mass flow rate, the various transport properties of the gases, and the 
geometry of the burner, and because these are largely invariant with pressure, 
these calculated heat-transfer coefficients are taken to be applicable at all 
pressures. This procedure, then, allows the solution of equation (C3) for T f 
at all other pressures. The resulting temperature curves are shown plotted in 
figure 8. Because of the extensive calculations and necessary approximations 
involved in obtaining these temperature curves, they must be considered approx- 
imate; nevertheless, they are believed to represent a fairly reasonable trend 
of flame temperature with increasing pressure. 
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APPENDIX D 


ESTIMATION OF THEORETICAL CO REACTION RATES 

Howard, Williams, and Fine (ref. 7) give for the global rate of oxidation 
of CO to form CO 2 , in the presence of water vapor 


5098 A Ccocifo^ (01 ) 

In a diffusion flame, combustion within the flame zone may reasonably be 
assumed to occur at an overall average equivalence ratio cj) = 1 ; hence, for 
the present CO-air flames, the initial composition can be taken to be 


xqq = 0.296 
xq 2 = 0.1 48 
x N2 = 0.556 


Hydrogen is almost always present in small concentrations because of virtually 
unavoidable contamination with water vapor, or by intentional addition. For 
most of the experimental runs, the CO cylinder with 99.3 ppm H 2 was used since 
no difference in results was observed between this concentration and the lower 
concentration of 1 0 ppm H 2 in the second cylinder. Therefore, for purposes of 
estimating the theoretical CO combustion rates, the concentration in the fuel 
gas was taken to be approximately 100 ppm or x H2 = 10~ 4 . In a stoichiometric 
mixture of CO fuel and air, this gives 

Xjj 2 = 2.96 X 10“ 5 


Taking x H20 ~ x H 2 ' an<3 substituting the given values of mole fractions into 

equation (Dl ) at T = Tf, the rate curves labeled r G in figure 10(a) are 
obtained. 

Several other expressions for the global rate of oxidation of CO also 
appear in the literature. (See, for instance, refs. 26 to 29.) For the same 

1 

flame conditions, these expressions predict rates from as much as 6— times the 

2 

r G values in figure 10(a) to as low as 1/5 these values, and some have somewhat 
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different temperature dependencies. The rate expression of Howard, Williams, 
and Fine (ref. 7) was selected here because it is more recent and because it is 
based on more extensive data. In addition, it predicts rates which are reason- 
able averages of those predicted by the other rate expressions. 

An alternate estimate of the rate of combustion of CO can be obtained from 
the elementary reaction steps 


CO 

+ OH + C0 2 + H 

(D2) 

CO 

+ O2 CO2 + 0 

(D3) 

CO 

+ 0 + M C0 2 + M 

(D4) 


The reaction rate constants k^ associated with these reactions are (ref. 8 ) 

k 2 = 4 x 1 0 1 2 e ~ 4026 / T 
k 3 = 1 .6 x 1 0 1 3 e -20634/T 
k 4 = 6 x 10 13 


with consistent units in moles, cm 3 , and sec. Reaction rates r 2 , r 3 , and r 4 
associated with reactions (D2) to (D4) were evaluated for both the 1 46-sccm CO 
flame and the 73-sccm CO flame as follows. The hydroxyl radical and oxygen-atom 
concentrations for use in reactions (D2) and (D4) were determined for equilib- 
rium conditions at T = Tf and for an equivalence ratio of hydrogen appropri- 
ate to the present system of 100 ppm H 2 in the CO fuel, <J> H = 1 0 “ 4 (i.e. r 

when the fuel is taken to be at an overall average $ of unity with respect to 
air). The results were then scaled for xq 2 =0.148 in the flame. The rate 

r 3 was calculated using the mole fractions for CO and O 2 given earlier. Com- 
parison of these calculated rates showed that « r 3 , and r 2 is roughly 
10 percent of r 3 for the 1 46-sccm CO flame, and roughly 2 percent of r 3 for 
the 73-sccm CO flame. Hence, a rough estimate of the theoretical reaction 
rates for CO can be written 

r T = r 2 + r 3 

where, for the 1 46-sccm CO flame, r T * 1.1 r 3 , and for the 73-sccm CO flame, 
r T * 1 . 02 r 3 . These rates were calculated at T = Tf and <J> = 1 and are 
plotted in figure 10(a). The close agreement between r T and tq for each 
of the two CO flow rates lends confidence that these rate predictions are of 
the proper order of magnitude. 
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APPENDIX E 


ESTIMATION OF DIFFUSION RATES 

To estimate the rate of diffusion of CO from the core of the flame into 
the reaction zone, Fick's law can be written as 


-► dc co 

J = -D 

dS- 


where 

J diffusional flux of CO 

D diffusion coefficient for CO (taken as CO into N 2 ) 

dc C0 

concentration gradient 

6l 

Expressing Cqq in terms of xqq by the perfect gas law, and writing the con- 
centration gradient in finite difference form, there results 


J = D 


( x CO)core “* ( x CO^edge 


RTf[ 


where the subscripts "core" and "edge” refer to the inner and outer edges of 
the flame zone, respectively, and L is the thickness of the flame zone (the 
same as the radiation path length in appendix C) . Taking (x co ) core = 1 and 
( x C0)edge = Of and writing the diffusional flux as a volumetric diffusional 
rate r D using the thickness of the flame zone to define the diffusional 
volume, there is obtained 


J pD 



Values for D were calculated using the method of Hirschf elder as described in 
reference 23, and L was taken from figure 18. 
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APPENDIX P 


ESTIMATION OF OXYGEN- ATOM CONCENTRATION 

In the combustion system of CO and air, the primary reactions involving 
oxygen atoms are (see, for instance, ref. 8) 


CO + o 2 -*■ co 2 + 0 

(Fl ) 

CO + 0 + M C0 2 + M 

(F2) 

0 2 + M 2 20 + M 

(F3) 

0 + N 2 -*■ NO + N 

(F4) 

N + 0 2 -*• NO + 0 

(F5) 


The reaction rate constants kj, in units of moles, cm 3 , and sec, for these 
reactions are given by 

k-j = 1.6 x TO 13 e -20634 /T 
k 2 = 6 x 1 0 1 3 


from reference 8, and by 

b = 3.56 x 10 18 - e -59386/T 
T 

k 3 r = 1 x l o 1 4 

k 4 = 1.36 x TO 14 e -37947 /T 

k 5 = 6.4 x 10 9 Te -31 45 / T 

from reference 9, where all rate constants are for the forward reaction as 
written, except the rate constant designated with a superscript r, which is 
for the reverse reaction. Since CO 2 and NO are stable reaction products which 
exit the flame zone, only the forward rates of reactions (Fl), (F2), (F4), 

and (F5) are considered; but since oxygen atoms are not a stable reaction 
product, the reverse rate of reaction (F3) is considered in addition to its 
forward rate. Employing a steady-state assumption for oxygen atoms in the 
flame, the rate of production of oxygen atoms is set equal to the rate of 
destruction of oxygen atoms. This equality can be expressed in terms of the 
rates of reactions (FI ) to (F5) as 
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APPENDIX F 


r-| + 2r3 + 1:5 = r2 + 2*3 r + *4 (F 6 ) 

where all rates r stand for the forward rates of the reactions as written, 
except the rate with the superscript r, which designates the reverse reaction. 
An order of magnitude analysis of equation (F6) for temperatures between 1600 K 
and 2200 K, for all pressures between 1 and 50 atm, and for mole fractions 
x CO = °« 296 , xq 2 = 0.148, x N2 = 0.556, with x M = 1.000, and with the assump- 
tion that atomic nitrogen is in equilibrium with molecular nitrogen, indicates 
that 

r 7 » 2r 3 
*1 » r 5 

r 2 = 0 ( 2 r 3 r ) for high oxygen-atom concentrations 
r 2 » r 4 

Hence, to a reasonable approximation at all temperatures and pressures in the 
flame, 

r-j « r 2 + 2r 3 r (F7) 

Substitution of the appropriate rate expressions into equation (F7) permits a 
solution for the steady-state concentration of oxygen atoms. In terms of mole 
fractions, this concentration is 


x O 


= -0.1 5x m + (o. O225x0 O + 6.565 - e - 20634 / T 


, 1/2 


X C0 X 0 2 


(F8) 


Taking xqq = 0.296 and xq 2 = 0.148, this simplifies to 


/ T V /2 

x O,ss = -0.0444 + h.971 x 1 0 -3 + 0.2876 - e -20634 / T I 


(F9) 


where the subscript ss is appended to designate that this oxygen-atom concen- 
tration is not an experimentally determined value, but is instead an estimated 
value derived from the present steady-state analysis. It should be pointed out 
that since the reverse of reaction (Fl ) was neglected in the analysis, the val- 
ues of x 0 ss obtained from equation (F9) probably represent an upper limit 
for x 0 which might be expected in the flame. The (F2) reverse reaction, 
which produces oxygen atoms, was also neglected; but because this reaction is 
highly improbable, its contribution to the concentration of oxygen atoms is 
insignificant. Values of x 0f ss evaluated from equation (F9) at the appropri- 
ate values of Tf and p, are plotted in figure 1 2 as *0, ss/^O^eq^- 
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3.06 i.d. x 0.05 wall 



Figure 1 . - Burner assembly and detail. All dimensions are in millimeters 
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(a) p = 1 atm. 

Figure 3.- Stability and extinction limits. 
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Figure 3.- Continued. 
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(a) 73-sccm CO. 

Figure 4.- Photographs of flame. 
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Figure 5.- Flame height as function of pressure. 
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Figure 7.- Molar emission index plotted against pressure for 
CO-air diffusion flames. 









CO flow rate, seem 

146 



Pressure, atm 


(b) Diffusional rate into flame zone 


Figure 10.- Concluded 
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Figure 14.- Required temperature for equilibrium molar emission index to 
equal experimentally observed molar emission index. 
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Figure 17.- Projected area of base of flame. 



Figure 18.- Flame thickness path length 
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